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Numerical studies of amorphous silicon show that the low- 
est 4% of vibrational modes are plane- wave like ( "propagons" ) 
and the highest 3% of modes are localized ("locons"). The 
rest are neither plane-wave like nor localized. We call them 
"diffusons." Since diffusons are by far the most numerous, 
we try to characterize them by calculating such properties 
as wavevector and polarization (which seem not to be use- 
ful), "phase quotient" (a measure of the change of vibrational 
phase on between first neighbor atoms), spatial polarization 
memory, and diffusivity. Localized states are characterized 
by finding decay lengths, inverse participation ratios, and co- 
ordination numbers of the atoms participating. 
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In region D, only diffusive propagation occurs over any 
meaningful distance, and the concepts of mean free path 
and wavevector lose usefulness. Although it may seem 
natural that the Ioffe-Regel crossover should be close to 
the mobility edge (Alexander, 1989), it is not true for the 
models we have studied. Indeed, Mott and Davis (1971) 
emphasize the non-coincidence. Both regions P and D lie 
in the extended (E) part of the spectrum. 

The non-coincidence of the E/L boundary and the P/D 
boundary means that the spectrum has three kinds of 
states. We find useful the terminology "propagon, diffu- 
son, and locon" given in Fig. ^|. The term "phonon" is 
avoided because of disagreement about what it means in 
a glass. 



I. INTRODUCTION 

Vibrational properties of disordered media have been 
reviewed by various authors, in particular Elliott and 
Leath (1975), Weaire and Taylor (1980), Visscher and 
Gubernaitis (1980), and Pohl (1998). Here we review and 
present new results in a program of numerical study of 
vibrations of amorphous Si. Among the new results not 
contained in earlier reviews are theoretical treatments of 
heat conductivity and thermalization rates in glasses. 

Harmonic normal modes of vibration can be rigorously 
classified as extended (E) or localized (L). In three dimen- 
sions the vibrational spectrum has sharp E/L boundaries 
("mobility edges") separating these two kinds of modes. 

There is another boundary, not sharp, which Mott 
called the "Ioffe-Regel limit" and which we call the "Ioffe- 
Regel crossover." This P/D boundary separates the spec- 
trum into a region with ballistic propagation (P) where 
wavevector is a reasonably good quantum number and 
a region with only diffusion (D) where wavevector can- 
not be defined but states are still extended. In region P, 
wave-packets can travel at sound velocity over distances 
of at least two or three interatomic spacings before scat- 
tering from disorder (Allen and Kellner, 1998). The dis- 
tance of ballistic propagation is the mean free-path t. 
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FIG. 1. Density of vibrational states from the 4096-atom 
model compared with data from Kamitakahara et al. (1987). 
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FIG. 2. Taxonomy of vibrations in glasses 

We instinctively seek additional labels, to replace the 
detailed classification scheme by wavevector and branch 
so useful in crystals. Here we attempt to characterize 
as completely as possible the dominant diffuson portion 
of the spectrum, but fail to find useful sub-categories. 
As already noticed by Kamitakahara et al. (1987), all 
diffuson modes of a given frequency have essentially in- 
distinguishable displacement patterns. 

We study amorphous silicon in harmonic approxima- 
tion, using a model we believe to be very realistic. As 
shown in Fig. [l], the E/L boundary is near the top of 
the spectrum, with only 3% of the modes localized. The 
P /D boundary is near the bottom of the spectrum, with 
only 4% of the modes ballistically propagating. Diffu- 
sons fill 93% of the spectrum. We do not think this is 
special either to our model or to amorphous Si. Other 
models for amorphous Si (Kamitakahara et al. 1987, Lee 
et al, 1991, Nakhmanson and Drabold, 1998) agree that 
the E/L boundary occurs near the top of the spectrum. 
Similar results are found for other glasses (Bouchard et 
al. 1988; Feldman and Kluge, 1995; Cobb, Drabold, and 
Cappelletti, 1996; Taraskin and Elliott, 1997; Carles et 
al, 1998), and model systems (Sheng and Zhou, 1991; 
Sheng, Zhou, and Zhang, 1994; Schirmacher, Diezemann, 
and Canter, 1998), with localized states occurring only 
near the top of the spectrum or in tails near gaps in the 
vibrational densities of states. For 3-d systems with ar- 
tificially large disorder, the mobility edge can be moved 
down to the middle of the spectrum (Canisius and van 
Hemmen, 1985; Fabian and Allen, 1996). The position 
of the P/D boundary near the bottom of the spectrum is 
widely accepted. 



II. LOW-FREQUENCY ANOMALIES 

Our study of normal modes by exact diagonalization 
on finite-size systems inherently lacks information at low 
frequency. For this reason, the present paper ignores the 
well-known but only partially understood low-frequency 
anomalies in vibrational properties of glasses. Our pre- 
vious work (Feldman, Allen, and Bickham, 1999) argues 
that the homogeneous network models we use probably 
contain no low-frequency anomalies. For completeness, 
we give here a brief catalog and point to sources for fur- 
ther information. 

"Two-level", or "tunneling" systems were introduced 
by Phillips (1972) and Anderson, Halperin, and Varma 
(1972) motivated by experimental discoveries by Zcllcr 
and Pohl (1971). The predictive strength of this concept 
is beyond question, but the physical objects the theory 
invokes remain elusive. A review was given by Phillips 
(1987). 

The phrase "boson peak" refers to a low frequency 
feature seen by Raman scattering (Stolen, 1970; Jackie, 



1981) in many glasses, which is correlated with the occur- 
rence of "excess modes" in specific heat and other spec- 
troscopies. There are many candidate explanations for 
these modes. One unified view, introduced by Karpov, 
Klinger, and Ignat'ev (1983) is called the "soft potential 
model." This model holds that glasses generically have 
anharmonic regions, modelled as double-well potentials. 
These give rise to two-level systems , relaxational be- 
havior, and quasi-localized or resonant harmonic normal 
modes. The last are a logical candidate for the excess 
modes. The subject was reviewed by Parshin (1994). 
Our models contain some quasi-localized modes at low 
frequencies, as is mentioned further in Sec. VIII. 

"Floppy modes" (Phillips, 1980; Thorpe, 1983) and 
the daughter concept of "rigid unit modes" (RUMs: 
Dove et al. 1996; Trachenko et al. 1998) refer to 
low-frequency modes which have zero restoring force in 
nearest-neighbor central-force models. Constraint count- 
ing algorithms provide methods of estimating the num- 
bers of such modes. They are expected to be quasi- 
localized in harmonic approximation, but intrinsically 
highly harmonic. Such modes probably do not play any 
important role in amorphous Si because of the overcon- 
strained coordination. 

Finally, experimental evidence shows that amorphous 
Si contains usually fewer two-level-type excitations than 
most glasses, and that samples with essentially zero such 
excitations can be prepared by treatment with hydrogen 
(Liu et al. 1997). 

III. THE MODEL 

Amorphous Si is an over-constrained network glass. By 
one usual definition, it is not a glass since it is not pro- 
duced by a glass transition upon cooling. Instead, thin 
films are condensed on cold substrates. When the film 
gets too thick, crystallization cannot be prevented. The 
absence of a glass transition can be attributed to good 
kinetic properties of a single-component system with one 
strongly-preferred bonding arrangement. To our mind, 
this just means that amorphous Si is simpler than most 
glasses, which for our purposes is more of an advantage 
than a disadvantage. With obvious caution required, we 
think that most of the properties we shall discuss can be 
regarded as typical of most glasses. 
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FIG. 3. Distribution of nth neighbors in the 4096-atom 
model of amorphous silicon. Fourth and fifth neighbors are 
shown in bolder lines. Note the small number of distant fourth 
or close fifth neighbors. The inset shows the distribution of 
bond angles averaged over all pairs of bonds with a common 
center and distances less than 2.6 A. 

For amorphous Si we use atomic coordinates gener- 
ated using the algorithm of Wooten, Weiner, and Weaire 
(1985). We have studied models containing 216, 1000, 
and 4096 atoms, contained in cubic boxes of side 16.5, 
28, and 44A, respectively, and continued periodically to 
infinity. All models are built using the Keating poten- 
tial (Keating, 1966) and then subsequently re-relaxed to 
a local minimum of the Stillinger- Weber potential (Still- 
inger and Weber, 1985). Different models differ in mi- 
nor details, both for ordinary statistical reasons and be- 
cause the algorithm was implemented slightly differently 
in each case. In this section we present structural prop- 
erties of a 4096-atom model. This model has larger dis- 
tortions from tetrahedral form than some of our other 
models, with 102 four-membered rings. Figure || shows 
neighbor distributions V n {uj) The distribution V\ for the 
first neighbor is calculated by measuring the distance to 
the closest neighbor of each atom; V2 is found from the 
distances to the second closest neighbor. Note that the 
first four neighbors distribute tightly within w 0.1 A to 
the crystalline first neighbor distance, 2.35A. The next 
12 neighbors in crystalline Si arc at 3.84A, a number fixed 
by the bond length (2.35A) and the tetrahedral bond an- 
gle (109.5°). In our model of amorphous Si, the closest of 
the next 12 (neighbor number 5) lies roughly between 3.2 
and 3.7A, while the farthest of these 12 (neighbor num- 
ber 16, not shown) lies roughly between 3.9 and 4.4A 
This reflects a flexibility in the bond angles, with values 
distributed between 90° and 125°, as shown in the inset 
to Fig. | The third set of neighbors in crystalline Si is 
12 atoms at 4.50A. This is determined by the fact that 
diamond structure has a dihedral angle of 60°, with all 



rings of the 6-member "chair" type. Rotation of the di- 
hedral angle to 0° ("boat" type rings) reduces the third 
neighbor distance to 3.92A Our model shows no gap 
at all between second shell (neighbors 5-16) and third 
shell (neighbors 17-28), consistent with random dihedral 
angles. 

The sum ^2 T^n over all n gives the radial distribution 
function g{r), plotted in Fig. |] and compared with the 
experiments of Kugler et al. The close agreement is one 
measure of the realism of the model. However, much of 
the structure of g(r) seems only to reflect atom density 
and nearest neighbor distance, so it may not be a very 
stringent test. 
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FIG. 4. Radial distribution function obtained by summing 
the nth neighbor distributions of Fig. ^ over all n, and divid- 
ing out a factor r 2 . The data are from a neutron diffraction 
experiment by Kugler et al. (1993). 



IV. VIBRATIONAL FREQUENCIES 

The other aspect of the model is the interatomic forces. 
We chose the model of Stillinger and Weber (1985) which 
is designed to work for the liquid as well as crystalline 
state. This required us to relax the coordinates from 
Wooten to a minimum of the Stillinger- Weber potential. 
The stability of the minimum is proven by the positiv- 
ity of all eigenvalues lu 2 of the dynamical matrix. The 
eigenfrequency distribution is shown in Fig. [l]. Qualita- 
tively satisfactory agreement is found with the neutron 
scattering data of Kamitakara et al. (1987). A simi- 
lar overestimate of vibrational frequencies is made when 
Stillinger- Weber forces are applied to crystalline Si, so 
we think the discrepancies should be attributed to the 
forces rather than the atomic coordinates. 
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V. LOCALIZED STATES 

The definition of a localized state is exponential decay 
of the eigenvector with distance from some center Rq: 



i{R n )\ oc exp(-|.R„ - #o|/&)- 



(1) 



This defines the localization length £j of the i-th normal 
mode, if the decay is observed. Fig. || shows selected 
modes, showing the very different character of modes 
from the D and the L portions of the spectrum. 



enough to locate the E/L boundary between 71 and 73 
meV; the precise location varies somewhat from model to 
model. We have made some experiments with artificially 
enhanced disorder, randomly scaling half of the masses 
by a factor of 5. This pushes the E/L boundary into the 
middle of the spectrum where it becomes more blurred 
by size effects. Localized states of "pure" amorphous 
Si are easily recognized. They are trapped in especially 
defective regions. This was discovered (Fabian, 1997b) 
by defining a local coordination number Z a , the number 
of neighbors of atom a at R a . The following arbitrary 
definition of neighbor suffices: 
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FIG. 5. Spatial decay of vibrational eigenvectors £j. Solid 
lines are guides to the eye. For each mode i the atom with 
the largest |e| 2 is located, and taken as the origin. The plot 
shows values of |e(7?„)| 2 averaged over spherical shells of \R n \ 
of width 8R = 0.2A. Results shown for w=50 meV, 71 meV, 
74 meV, and 87 meV, are averages over 3, 4, 2, and 2 modes 
respectively. Modes at 50 meV and 71 meV have mean square 
eigenvector near l/(iV=4096) throughout the cubic cell of 
length 44 A, with enhanced values only in spatial regions of 
small measure. These modes are extended, as are all the 
modes below the mobility edge at 73 meV. Modes at 74 meV 
and 87 meV have mean values of \e(R n )\ 2 « exp(-2|.R„|/£) 
with localization length £ = 4.7A and 1.5A respectively. 

There is still controversy concerning the location of 
the mobility edge in glasses. Unfortunately, most exper- 
iments shed little light, since measured spectral proper- 
ties, being averages over a macroscopic region, do not dif- 
ferentiate between localized and delocalized states. Heat 
conductivity k(T) is the property most strongly affected 
by localization. We think the measured k(T) strongly 
supports our placement of the P/D and E/L boundaries 
near the lower and upper edges of the spectrum. This 
will be discussed in the next section. 

If eigenvectors are calculated, then no special tricks are 
needed to locate the theoretical mobility edge in models 
of amorphous silicon. A model with 216 atoms was large 



where w(r) is 1 for r < 2.35A, for r > 3.84A, and 
continuous and linear in between. This gives an aver- 
age coordination of 4.7 neighbors. The "mode average 
coordination number" is defined as 



Zj = Z a \ti{a)Y 



(3) 



Most modes have Zi near average, but localized modes 
mostly collect at regions with significantly higher coordi- 
nation, as shown on Fig. ^. 
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FIG. 6. The mode average coordination number (Eqn. |^), 
and the participation ratio (on a logarithmic scale) , calculated 
for a model with 1000 atoms. 

We have seen the mobility edge in many independent 
calculations. 

1. The participation ratio pi (Bell and Dean, 1970), 
plotted on Fig. ^, gives the number of atoms on 
which the mode has significant amplitude. For the 
E=P+D part of the spectrum, the value hovers 
near 500 for a model with 1000 atoms. 
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2. The diffusivity Di drops to zero at the mobility 
edge. This is shown in the next section. 

3. Sensitivity to boundary conditions is also discussed 
in the next section. 

4. Level spacing distributions have been computed for 
diffusons and locons (Fabian, 1997b). As expected, 
diffusons obey Wigner-Dyson statistics, while lo- 
cons obey Poisson statistics. 

5. Eigenvector self-correlation functions were carefully 
studied in a 4096-atom model by Feldman et al. 
(1998). 

6. Many quantities qi which can be evaluated for each 
mode i seem to depend only on Wi for diffusons, but 
become mode-specific for locons. The ones we have 
looked at are: 

(a.) The mode-average coordination number, plot- 
ted on Fig. ||. 

(b.) The phase-quotient parameter, discussed in 
section VII. C. 

(c.) Bond-stretching parameters (Fabian and 
Allen, 1997). 

(d.) Griineisen parameters for both volume and 
shear deformations (Fabian and Allen, 1997 
and 1999). 



latter model is only justified if the mean free-path £q — 
vqtq is longer than a wavelength. 

It was noticed by Birch and Clark (1940), and by Kit- 
tel (1948) that in glasses k(T) at T >20K could be in- 
terpreted as the specific heat C(T)/V multiplied by a 
temperature-independent diffusivity D of order a 2 cj_o/3 
where a is an interatomic distance. In the phonon-gas 
model, this would correspond to £ « a, too small to jus- 
tify use of the model. The success of this observation 
implies that the dominant normal modes in a glass are of 
the D variety, not P because P implies £ ^> a, and not L 
because L implies D — until anharmonic corrections are 
added which make D depend on T. This successful (and 
we believe, essentially correct) interpretation lost favor 
after Anderson localization was understood, because a 
misconception arose that the P/D boundary (which cer- 
tainly lies low in the spectrum of a glass) should lie close 
to the E/L boundary. 

Our numerical calculations of Di are shown in Fig. [?]. 
Also shown are values of Di from a formula of Edwards 
and Thouless (1972), 



Di = L 2 Alo, 



(0) 



where Awi is the sensitivity of the eigenenergy to a twist 
of the boundary condition. We have simply used for Au>i 
the change in Ui when boundary conditions are changed 
from periodic to antiperiodic. The actual definition is 
probably 



VI. HEAT CONDUCTIVITY AND DIFFUSIVITY 

Theoretical interpretation of heat conduction of glasses 
has been contentious. Below the "plateau" region (typ- 
ically 5-30K) it is agreed that heat is carried by ballis- 
tically propagating low-frequency modes (the P region 
of the spectrum). Above the plateau, k(T) rises, ap- 
proaching at room temperature a constant value which is 
typically smaller than the crystalline value (a decreasing 
function of T at room temperature.) A rigorous conse- 
quence (Allen and Feldman, 1993) of the Kubo formula 
(Kubo, 1957) and the harmonic approximation is the re- 
lation 



«(T) 



^2 C(hLJ i /2k B T)D i , 



(4) 



where C(x) is the specific heat of a harmonic oscillator 
(x/ sinh(a;)) 2 and Di is the "diffusivity" of the i-th. nor- 
mal mode of frequency u>i, given by 



Di 



ttV 2 v 



\Sij\ 2 S(uJi - Uj), 



(5) 



where S%j —< i\S\j > is the intermode matrix element of 
the heat current operator. Eq. || also emerges, with Dq 
equal to VqTq/3, from the Peierls-Boltzmann phonon- 
gas model (Gurevich, 1986) of transport in crystals. The 



Aw, 



lim 

0^0 



(7) 



where Auji (0) is the shift when the phase is twisted by <fi. 
Antiperiodic boundary conditions correspond to <j) = 7r, 
while <j) = 2n returns to periodic boundary conditions 
with Aa>i(27r) = 0. Therefore our calculation, which is 
the only one easily accessible for us, gives a probable 
upper bound to Di for each mode i. Inspection of Fig. |?j 
shows that with this interpretation, the two calculations 
agree reasonably well. Both go to zero at the mobility 
edge, and both become large and ragged in the P region 
below 10 meV. The raggedness comes from the sparseness 
of the eigenstates at low oj, and the large values reflect 
the onset of ballistic propagation. In the D region above 
12 meV, values have collapsed to the range of 1 mm 2 /s, 
which corresponds to W£>a 2 /3, with u>d set to 50 meV 
and a = 2A. This diffusivity is well below any value that 
could be allowed in a phonon-gas picture, and agrees with 
the measured k(T) (Allen and Feldman, 1990; Feldman 
et al, 1993). The peak of Di around 33 meV corresponds 
to a feature in the "phase-quotient" that will be discussed 
in section VII. C. Similar results for vitreous Si02 have 
been reported by Feldman and Kluge (1995). 
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FIG. 7. Diffusivity versus frequency (Feldman et al., 1993, 
Fabian, 1997b) calculated for a 1000 atom model by Eq. [s| 
and for a 4096 atom model by Eq. ^| Mobility edges, shown 
as vertical lines, are slightly different in these two models. 

When k(T) is calculated from Eq. ||, using the values of 
Di from Fig. [?], the results, shown in Fig. || agree roughly 
with the data at higher temperatures. At low tempera- 
ture it is necessary to have an additional source of heat 
current, the ballistically propagating long- wavelength 
modes. In Fig. |[ this has been added in a thoroughly ad 
hoc fashion. We have simply assumed a Debye spectrum 
for the modes with uj < luq ~3 meV, and a tempera- 
ture independent diffusivity D(u>) = Dq{ojq/ijj) 2 . There 
is no theoretical justification for this. In principle, the 
temperature-independent diffusivity caused by glassy dis- 
order should take the Rayleigh (luq/uj) 4 form at low u>, 
and one needs a stronger type of scattering, inelastic, and 
therefore T-dependent, to match the data. However, the 
(loq/u)) 2 behavior has been seen at intermediate frequen- 
cies, both experimentally (Sette et al. 1998) and numeri- 
cally (Dell'Anna et al. 1998; Feldman et al. 1999), so we 
have used this simpler fitting proceedure to match the 
data. 
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FIG. 8. Thermal conductivity k, versus temperature T mea- 
sured for amorphous silicon (circles: Pompe and Hegenbarth, 
1988; squares: Cahill et al., 1989; diamonds: Cahill et al., 
1994). The bold line is a calculation from the diffusivity val- 
ues shown in Fig. FJ for a 1000 atom model using Eq. m The 
thin dashed line is a calculation from modes at low energy 
assuming that their diffusivity obeys an co~ 2 law. The bold 
dashed line is the sum of these two contributions. 

Our most important conclusion is that the "re- 
increase" of thermal conductivity above the plateau re- 
gion is attributable to heat carried by "diffuson" modes in 
much the way imagined by Birch, Clark, and Kittel, and 
that the plateau is a simple crossover region, not requir- 
ing any new physics to explain. In particular, we believe 
that "excess modes" (also known as a "Boson peak" ) is 
not a necessary ingredient to explain the plateau. Amor- 
phous silicon seems to lack these "excess modes" but still 
to have a plateau. 



VII. THERMAL EQUILIBRATION 

There is some evidence suggesting that vibrations in 
glasses, if disturbed from equilibrium, may return very 
slowly. For amorphous silicon, experiments were reported 
by Scholten and Dijkhuis (1996) and by Scholten, Aki- 
mov, and Dijkhuis (1996). Our investigations show that 
if the disturbance is not too large and is purely vibra- 
tional, then the rate of return to equilibrium should be 
as fast, if not faster, than in a corresponding crystal. Sur- 
prisingly, we find that this is true both in the locon and in 
the diffuson portion of the spectrum, contradicting a view 
(Orbach, 1996) supported by fracton theory (Alexander, 
1989) that localized vibrations must equilibrate slowly. 
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FIG. 9. Decay rates calculated by anharmonic perturba- 
tion theory for amorphous Si. The inset shows the process 
diagrammatically. 

In thermal equilibrium, the harmonic vibrational 
eigenstates have average population given by the Bosc- 
Einstein distribution. Fabian and Allen (1996) used an- 
harmonic perturbation theory to calculate the inverse 
lifetime or equilibration rate 1/r by which a vibrational 
state returns to the Bose-Einstein distribution if driven 
out of equilibrium. Their results are shown in Fig. [|. 
The validity of the perturbation theory is confirmed both 
by the smallness of the ratio 1 / lot and by a direct test in 
the classical regime using molecular dynamics by Bick- 
ham and Feldman (1998) and Bickham (1999) 

It can be seen in Fig. ^ that no change occurs in the 
size of 1/r at the mobility edge. A careful treatment of 
locons does not support the idea that they thermalize 
slowly. Their anharmonic thermalization rates (Fabian, 
1997a) are as fast, or even faster, than diffusons, and 
also comparable to, or faster than, corresponding ther- 
malization rates of vibrations in crystals. The source of 
the misconception of slow equilibration is the idea that 
direct decay of a locon into two locons should be negligi- 
ble. This idea fails because, unlike for example, band-tail 
electronic states in amorphous Si, the vibrational states 
are not at all dilute. Slow thermalization rates (forbidden 
by all theories we understand) could be tested by looking 
for a contribution from thermal vibrations to attenuation 
of very high-frequency sound (Fabian and Allen, 1999). 
If the thermalization rate is extremely slow, this contri- 
bution to the attenuation would be greatly enhanced. 



VIII. RESONANT MODES 

Inspection of Fig. ^ shows that a few modes in the D 
region and somewhat more in the P region have anoma- 
lously small participation ratios, of order 100 out of the 
1000 atoms available. These states are not exponen- 
tially localized (Fabian, 1997b; Feldman et al, 1999) 
but are temporarily trapped in regions of peculiar coor- 
dination, from which they can tunnel into the continuum 
of extended states. Such states were first reported by 
Biswas et al. (1988) in a small model similar to ours; a 
model with larger numbers of 3- and 5-fold coordinated 
atoms had more such modes, which were speculated to 
bear some relation to the "floppy modes" of Phillips and 
Thorpe. Such modes were studied in detail by Schober 
and coworkers (1988, 1991) and Oligschleger and Schober 
(1999). 

We have recently argued (Feldman et al, 1999) that, in 
our (mostly 4-fold coordinated) models of amorphous Si, 
such states tend to disappear as the size of the model gets 
bigger, presumably because each such mode is trapped 
only in a very specific peculiar region. As the number 
of atoms in the model increases, so does the number of 
peculiar regions, but if each resonant mode is trapped in 
only one region, the fraction of time spent outside that 
region increases because the volume outside that region 
has increased. On the other hand, such modes, especially 
the ones in the P region, may be more pronounced in real 
amorphous Si and other real glasses than they are in the 
models we study. This is because our models are spatially 
homogeneous on scales greater than 4 A, while real glasses 
may have mesoscopic defects such as voids which would 
attract more such modes. 

Fabian and Allen (1997, 1999) found that the reso- 
nant modes have giant (« —40) Griineisen parameters 
7j. These parameters measure the sensitivity of loi to 
macroscopic strain. In a glass (just as in a crystal where 
positions of atoms are not all fixed by crystallographic 
symmetry) strains cause not just a homogeneous shift 
of atomic coordinates, but also an additional local re- 
laxation, which turns out to be particularly large in just 
those peculiar regions where the resonant modes are tem- 
porarily trapped. Anomalously large values of 7, play an 
important role in explaining the anomalously large and 
sample-dependent measured low-T thermal expansion of 
glasses, and also should show up in enhanced contribu- 
tions to the attenuation of high to sound waves at higher 
T. 



IX. CHARACTER OF DIFFUSONS 

The most important property which distinguishes dif- 
fusons is their intrinsic diffusivity Di with values of order 
ujdo 2 /3. If wave-packets were constructed in such a way 
as to be approximately monochromatic, and simultane- 
ously localized at the center of a cell with a reasonably 
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small radius (perhaps 6-8 A) , then we believe that no mat- 
ter how well directed such a pulse was at t = 0, the cen- 
ter of the pulse would hardly move, and the radius would 
evolve as < r 2 >= 6Dt for all times until reaching the 
cell boundary where it would interfere with its periodic 
image. Unfortunately, a 44A cell is only marginally big 
enough, and computational difficulties have so far pre- 
vented us performing this experiment. Here we describe 
our efforts to find other ways to characterize diffusons. 



A. wavevector 

At the P/D boundary, wavevectors become ill-defined. 
Fig. [Ic] shows a test. The squared Fourier weight is 
defined as 



Wi(Q) 



,iQ-R(t); 



(8) 



where Q is chosen as (2Tr/L)(h, k, I) with integer h, k, I so 
that the periodic images interfere constructively. We de- 
fine Wi(Q) as Wi(Q) averaged over spherical shells of wave 
vector of width 0.2 x 2n/L. The value Q = 9.2 x 2n/L 
corresponds to neighboring atoms being completely out 
of phase. The 51 meV diffusons show a peak Fourier con- 
tent near 8x2ir/L, but the peak height is less than twice a 
"background" value found at larger Q which dominates 
the behavior. There is no ballistic character to these 
modes. 



B. polarization 

Diffusons have no wave vector. Not surprisingly, 
they also lack a polarization as is shown in Fig. \VU. 
Propagons, by definition, have a wavevector. The na- 
ture of the propagons in the 4096 atom model was ex- 
amined by Feldman et al. (1999). As shown there, the 
modes near u> =3.5 meV have well-defined TA character, 
with the smallest possible wave vector Q = 2ir/L. Fig 
|TT| shows that these modes have only limited preference 
for a direction of polarization. Similarly, the mode at 
7.2 meV has Q = 2ir/L and LA character, but not much 
polarization. 
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FIG. 11. Area-preserving projection normalized unit polar- 
izations ii(R(£)) onto circles (one circle for each hemisphere.) 
Three modes i are shown, with each of the 4096 components 
showing as a dot. 
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FIG. 10. Fourier weights calculated from Eq. ^ for diffu- 
sons at 50 meV and 71 meV, and locons at 74 meV, averaged 
over 3, 4, and 2 modes respectively. 
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FIG. 12. The spatial fall-off of e(fi) ■ e(fj) weighted by nj 
and averaged for the same 3 modes at 50meV and 4 modes at 
71meV as in fig. 
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Polarization directions of diffusons wander uniformly 
over the unit sphere, and one may ask what is the spa- 
tial range of decay of polarization memory in a mode. 
This is shown in Fig. ^[ In a crystal, one has sur- 
faces in Q-space where modes have constant frequency 
lu(Q) =const. A single eigenstate is some arbitrary lin- 
ear combination of the degenerate Bloch waves on this 
surface. If the surfaces were spherical and the linear 
combinations random, we would expect that e(fj) ■ 
would fall off spatially as cos {Qrij)/rij. Fig. [12] shows 
that for diffusons at 50 meV, some polarization memory, 
but much less than expected for a crystal, persists out to 
12 A, while higher to diffusons lose polarization memory 
more rapidly. 

C. phase quotient 




10 20 30 40 50 60 70 80 90 
FREQUENCY (meV) 



FIG. 13. Phase quotient versus energy for a 1000 atom 
model. 



The "phase quotient" fa was defined by Bell and 
Hibbins-Butler (1975) as 



E< a ,b> 13(a) I' 



(9) 



and is plotted in Fig. [13| For low u>i, values near 1 indi- 
cate that nearest neighbor atoms (the only ones summed 
in Eqn. [)]) vibrate mostly in-phase, while for high u>i, 
values near -1 indicate that nearest neighbors vibrate 
mostly out-of-phase. Like so many other properties, this 
depends only on uji until the E/L boundary is crossed, 
but is very mode-specific for locons. The sharp rise at 
LUi w 29 meV is interesting, and may help explain why 
at the same frequency Di (Fig. [?]) has a sudden rise. 
In crystalline Si, at approximately the same point in 
the spectrum, the TA branch ends and the density of 
states has a local minimum. Thus 29 meV marks a point 
in the spectrum where diffusons change character from 
bond-bending (somewhat TA-like) with relatively high 
frequency because of large phase difference from atom 
to atom, to bond-stretching (somewhat LA-like) with 
not such a large phase difference but an equally high 
frequency because the bond-stretching forces are bigger 
than the bond-bending forces. Apparently the latter kind 
of mode has greater diffusivity by a factor 2 or more. In 
a crystal we attribute this to a larger group velocity of 
the LA branch and a smaller density of states for de- 
cay by elastic scattering. Neither of these properties can 
be properly invoked for diffuson modes in a glass, but 
apparently similar physics is somehow at work. 



X. SUMMARY 

Since 95% of the states in amorphous Si, and proba- 
bly many other glasses as well, are diffusons, we should 
understand their properties. All diffusons at a given fre- 
quency u> seem essentially identical. As u> changes, their 
properties evolve, mostly smoothly, but the sudden jump 
in diffusivity and in phase quotient at 29 meV shows that 
not all variability is lost. 
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